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We suggest that extension of the Racah seniority description of strongly interacting fermions in the 
nuclear shell model is directly generalizable to describe pairing of atoms in cold Fermi systems. We 
illustrate this by the fermionic pairing in the much studied cold two-component gas of 6 Li atoms. Our 
pairing interaction is two orders of magnitude stronger than that used in the usual BCS approach. 
We also explain why the Racah scheme is less applicable to nuclei, and discuss the similarities of the 
strongly-coupled matter in cold fermion systems and the new form of matter found in RHIC close 
to T c . 



I. INTRODUCTION 

Ultracold atomic gases have become a medium to re- 
alize novel phenomena in condensed matter physics and 
test many-body theories in new regimes [1]. Recently 
Regal et al. [2] observed a preexisting condensation of 
fermionic atom pairs in the region of magnetic field B 
above the Feshbach resonance, generally referred to as 
the BCS side. Whereas Regal et al. used a trapped gas 
of fermionic K atoms, Zwicrlein et al. [1] found with Li 
atoms a much larger condensation of fermionic atoms in 
the same BCS region. Such a large condensation of es- 
sentially zero momentum molecules on the BCS side is 
difficult to understand in the weak coupling BCS theory 
following from the Fermi-Thomas approximation usually 
employed in theoretical discussions. 

By staying in the harmonic oscillator representation 
(rather then approximating by plane waves), we uti- 
lize the large degeneracy of atoms within a major shell 
N. This gives a factor of about 60 increase in pairing 
strength. Pairing is between time reversed states m and 
-m, giving a simple explanation for the zero total mo- 
mentum of each molecule in the condensate. 



II. GENERALIZED SENIORITY 

Racah [3] showed that for a <5-function potential the 
ground state solution for the J™ shell was given by cou- 
pling all pairs to J = 0. This is reasonable, because with 
fermions, antisymmetry requires that 
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so that fermions in one pair will have zero average in- 
teraction with those in another pair. Thus, the energy 
of a J™ configuration is the sum of pairing energies (the 
seniority quantum number being the number of unpaired 
particles) . 

Even though the seniority representation gave the ex- 
act answer only for an attractive zero range interaction 
in the J™ configuration, it could easily be extended to 
other degenerate configurations [4] . With a pure pairing 
interaction 
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with matrix elements assumed equal for all angular mo- 
menta I in a major shell of a principal quantum number 
N, the collective wave function 



* = 



1 

V / 
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where the states are the paired states of azimuthal quan- 
tum numbers m and -m and Q is the number of pairs 
(m, —tn), is moved downward in energy with full trace of 
the secular matrix. This is just the bosonization of the 
fermion pairs. We have chosen an example of the N = 15 
shell, with a trap energy 15 tkjQ, of a cold Fermi gas with 
6 Li atoms in mind, with all levels through N = 15 being 
filled. The number of atoms is 1632. 

We first make an estimate of the degeneracy. Since in 
the harmonic trap 
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then in dimensional analysis, the Slater integral 
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and the total number of states [5] in shell N is (N + 3/2) 2 
in our assumption where all Slater integrals, both for 
1 = 1' and I ^ V contribute equally, we find a pairing 
energy* 
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"The coupling constant in eq. (2) is given in terms of the 
negative scattering length a by G = 47rft 2 |a|/m. 
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so that 



= 4tt N 



3\ a 



= 51 — 
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where a„ 



(ft/m^) 1 / 2 



We compare this with Heiselberg's eq. (16) [6], 
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with the last shell Nf — 15. 

A shell-model calculation, without use of the Thomas- 
Fermi approximation, in the Nf = 12 shell gave 
1.14|o|/a osc to compare with Heiselberg's 1.06, showing 
the Thomas-Fermi approximation to be good. 

Because of the large factor between our rough esti- 
mate, eq. (7) and Heiselberg's, we undertook a numeri- 
cal diagonilization of the pairing interaction (eq. (2)) in 
the TV = 15 shell without assumption of equal matrix 
elements using harmonic oscillator wave functions. The 
result was that the factor 51 in eq. (7) was replaced by 
58.6.^ We note that this latter number is 4n times 58% 
of the trace of the matrix we diagonalized. Using the 
Brown-Bolstcrli model [7], in which off-diagonal matrix 
elements are approximated by 



Mij = y/Ma Mj 



(9) 



the factor would be An times the trace of the matrix, so 
that, within factor of about 2, the Brown-Bolsterli fac- 
torisable model, which involves calculation of only the 
diagonal matrix elements, can be used. This should fa- 
cilitate calculations involving many shells. 

From the similarity of the mathematics with the 
Brown-Bolsterli model, we see that it is the degeneracy 
of the configurations within the shell of a given princi- 
ple quantum number which gives us our factor of about 
60 as compared with weak coupling BCS. In fact, as can 
be seen, the usual BCS treatment gives a rather weak 
pairing, with a gap 
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aside from a factor of about unity. In general, kf is taken 
as y/2Nf /a osc , where Nf = 15 in our case. This pair- 
ing is a surface effect, holes being made over a region 
of momenta about the Fermi surface so that particles 
of equal magnitude, but opposite momenta, can scatter 
with them. 

In our case of a volume effect, the gap begins linearly 
with a, increasing rapidly with increasing degeneracy. 



^This is an underestimate because backward-going ladders 
(ground-state correlations) which we have not included give 
additional attraction. 



With the large binding energy of eq. (7) we can eas- 
ily increase A so that it cancels the trap energy for the 
Nf — 15 shell. In fact, we can make the pair of fermions 
into a boson, following a scenario somewhat similar to 
that suggested by Falco and Stoof [8]. The molecular 
binding energy of the boson is 
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We can provide the total binding energy from the mean 
field energy on the right hand side of eq. (7) : 
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which gives 



= 1.75. 



Using the kf = y / 2Nf/a osc [6] this gives 
' =0.32. 
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For 8x 1632 atoms (kf |tz| ) 1 would be about half of this. 

Our factor 60 increase is achieved chiefly by correla- 
tion in angle. By pairing the m and -m projections of 
angular momentum and summing over m one gets using 
the addition theorem for angular momenta 



*L=0 



V2T+T Pi (COS 0i 2 ) 



for the angular dependence of the paired state [4] . The 
factors \/2l + 1 were used in the Brown-Bolsterli model 
[7] which was formulated in the representation of good 
angular momentum; in the m representation the matrix 
in eq. (2) with all matrix elements equal is equivalent. 

Of course, the S- function interaction sets Pi (cos #12) = 
1 so that the factor (2/+ 1) from squaring 5'l=o resulted 
in Racah's seniority scheme. We gain additionally by 
including the off-diagonal Slater integrals. 

In fact, our interaction is increased by the two-body 
correlations, so that kf, which is determined by the to- 
tal number of particles divided by the volume of the 
trap; i.e., by the average density, should still be given 
by y/2Nf /a osc . We can, therefore, rewrite our eq. (7) for 
A as 



A 



= n kf \a\ 
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neglecting the 3/2 as compared with N, so that the N 
dependence is taken up in kf. In this way we do not 
have to know a osc . The a/2-/V/ times 11 gives the factor 
of about 60 in our enhancement in the N — 15 shell. 

It should be noted that for the spherical trap we 
have assumed that 34 out of this 60 comes from the 31 
Z = 15 pairs and that our seniority solution is exact 
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for a S- function interaction [3]. In the experiments the 
traps are not spherical, typically being cigar shaped with 
quite large ratios of frequencies. In nuclear physics, al- 
though with lesser ratios, one uses the Nilsson model, 

with uj = (uj x uj y LOz) 1 ^ . Since, (r 2 ) 1 ^ , in cq. (4), goes 

as oj 1 I 2 this means that the dependence of (r 2 ) 1 ^ on 

any particular Wj, say lu z , goes as luz 1 ^ 6 . Now in a cigar- 
shaped trap there will still be the correlation in 612 en- 
hancing the densities of time-reversed orbits m and -m, 
but, especially at the equator, the distance r from the 
center of the trap to the edge of the cigar will be smaller 
than that of a spherical trap with the same number of 
m values. Thus, the distance r #12 will be shorter than 
in the spherical trap, so the enhancement in probability 
will be greater. However, this should not be a large ef- 

1 /9 

feet because of the weak dependence of (r 2 ) N[ on any 
particular Wj. 

For small values of \a\, the particles in the levels be- 
low N — 15 will interact with those in JV = 15 so as 
to push the level up (although in the middle levels they 
will push the levels above up and the those below down) . 
However, the attraction goes as VN, as can be seen from 
eq. (7), whereas the trap energy goes linearly with N, 
which decreases faster than the attraction with decreas- 
ing N, so that the lower levels will have already gone 
into molecules, the atoms from shell Nf being the last to 
make the transition. 

We have not done dynamics, but we believe that for 
very small \a\ there will be a collisionless regime, in which 
the RPA vibrations of Bruun and Mottelson [9] occur. 
The binding energy of the system, the sum of bubbles 
which make up the RPA vibrations or zero-point energy 
of the vibrations [10] can be incorporated, as noted by 
Bruun and Mottelson, as a mean field which slightly in- 
creases the oscillator energy spacings. But then as a be- 
comes larger in magnitude the inner shells start losing 
their support as the attractive field from the interaction 
reaches the trap energy in magnitude and the shells start 
collapsing in what becomes a hydrodynamical (strongly 
interacting) regime. The end result, already early on for 
negative values of (kja)^ 1 is that the strongly bound 
molecules are formed. The higher shells of filled atoms 
may remain, but their hu>o will be lowered by the attrac- 
tive mean field. 



III. SENIORITY IN NUCLEI 

We explain why nuclei, for which Racah introduced 
seniority, do not look like the cold Fermi gases we have 
been discussing. Although seniority was very useful in 
discussion of the energies of nuclei, it has turned out 
that the short range attraction in nuclei is totally oblit- 
erated by the short-range repulsion, so that the effective 
interaction in s-wave is only the long-range part of the po- 
tential, from about 1 fm outwards. This is most clearly 



understood in the Moszkowski-Scott separation method 
in the calculations of Holt and Brown [11]. In the calcula- 
tion of nuclear interactions, as done by Kuo and Brown 
[12] which used this separation method, the inner part 
of the s-wave two-body interaction up to separation dis- 
tance d ~ 1 fm is simply removed, leaving the long-range 
part of the two-body interaction Vi to be used in shell- 
model calculations. Holt and Brown showed that this 
Vi is the configuration space form of Vi ow -k, which is the 
more effective effective field theory (MEEFT) interaction 
and contains all experimental information in the nucleon- 
nucleon scattering phase shifts [13]. This long range Vi 
is clearly more suitable for shell-model calculations than 
for pairing. The latter in nuclei is of the weak coupling 
BCS type and takes place only in the last shells, near the 
Fermi energy. 

IV. RELATION TO RHIC PHYSICS 

Edward Shuryak [14] has often emphasized that the el- 
liptic flow seen by O'Hara et al. [15] of the cold Fermi liq- 
uid, when released from the trap, is similar to that formed 
at RHIC, although there are many orders of magnitude 
difference in them (e.g microseconds versus fermi/c). Our 
description of the pairing involves the same wave func- 
tion * (eq. (3) as that used by Brown, Lee, Rho and 
Shuryak [16], (Appendix A. in [17]) , although at RHIC 
the bosonization of quarks and antiquarks into chirally 
restored mesons above T c is involved, rather than fcrmion 
pairs. There the quark and antiquark masses from lat- 
tice gauge calculations are found to be ~ 1 GeV, and 
the 7r and a meson masses must be brought to zero at 
T = T c , in the chiral limit, which requires extremely 
strong interactions. The attractive interactions are built 
up enormously through the degeneracies. 

Thus, although the Brown-Bolstcrli model for giant 
resonances in nuclei suggested the use of the degenera- 
cies both here and in RHIC, they are immensely strong, 
because the numbers are larger and the degeneracy is not 
partially broken by spin-orbit coupling, as it is in nuclear 
physics. Inclusion of the degeneracy brings about the 
strong coupling. 

V. DISCUSSION 

Our result is consistent with the results of Carlson et 
al. [18] who found A to be about half of the Fermi en- 
ergy in quantum Monte Carlo studies of superfluid Fermi 
gases. In [17] Brown et al. discuss fluid formed in RHIC 
in terms of Nambu Jona-Lasinio theory. They find that 
the interactions above T c are much stronger, because of 
the large degeneracy, than below T c , where the unper- 
turbed spectrum of quarks and quark-holes is continuous 
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as in the Fermi gas.* 

We believe the experiments of Regal et al. [2] and 
Zwierlein et al. [1] found pairs of atoms with properties 
such as our scenario would produce. First of all, our 
molecules are formed on the BCS side so one would ex- 
pect large condensate fractions there when the tempera- 
ture was lowered. They find a large fraction of molecules 
to be in a zero-momentum state after fast ramping in 
magnetic field over to the BEC side, "this means that 
the nearest neighbors had opposite momenta" [1]. In 
our seniority scenario, atoms of time-reversed states m 
and — m are paired, states of equal in magnitude and op- 
posite in direction angular momenta. In the long-range 
Cooper pairs discussed in other works one would expect 
the transfer into a tightly bound molecular state would 
randomly pick out one of the nearest neighbors, result- 
ing in a thermal molecule. Zwierlein et al. [1] regard their 
high condensate fraction as evidence of condensed atomic 
pairs in the BCS sector which are smaller in size than the 
inter atomic distance and, therefore, molecular. The size 
of our pairing modes 'J, eq. (3), decreases roughly in- 
versely with increasing f2, the number coupled together, 
forming a highly compact molecular state. 

The condensate drains particles from the Fermi sea on 
the BCS side as we have outlined, but only after they are 
tightly paired. They go into zero-momentum molecules 
so that they are easily Bose condensed. As noted in [1] 
both there and in [20] and in [21] a reduction in the size of 
the atomic cloud was observed as the Feshbach resonance 
was approached from the BCS side. 

We believe that our generalized seniority scenario ex- 
plains the important features of the qualitative change in 
pairing phenomena demanded by the Regal et al. [2] and 
Zwierlein et al. [1] results. 

Our conclusion is that the generalized seniority de- 
scription gives a new strong-coupling pairing which de- 
scribes important features found in experiments. 

We end by noting that Leggett's ansatz [22] is 



with the quasi molecular wave function as a variational 
parameter. The </>, common to all pairs, is 



<f>(ri <J 1 ,r 2 a 2 ) = -j= (T I 



IT) <P(\ri-r 2 \ 



It reduces to a Bose condensation of tightly bound di- 
atomic molecules in his limit £ = (fcja) -1 — > oo and to 
the standard BCS result in the limit £ — > — oo. This <f) is 



"'"In a longer paper we shall show schematizing further the 
schematic model of a Cooper pair [19] that it is possible to go 
smoothly from the weak coupling BCS to our strong coupling 
one, as the degeneracy is introduced. 



well mimicked in the harmonic oscillator representation 
by our ^ of eq. (3), which is made up out of the trap har- 
monic oscillator wave functions and, therefore, preserves 
the degeneracy in a major shell. 
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